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We show that both hard perturbative and soft non-perturbative QCD effects are 
important in hadroproduction of charmonium. Observed xp and distributions 
of J/tj) in fixed target experiments can be well described by our Monte Carlo model 
based on perturbative QCD and the Color Evaporation Model. 

Aiming at an understanding of non-perturbative QCD (non-pQCD), we 
study effects of soft color interactions on charmonium production. This is a 
particularly interesting case, since the charm quark mass rric is a large enough 
scale for perturbative QCD (pQCD) to be used as a basis, but small enough 
for non-perturbative effects to be of importance. This provides an interesting 
interplay between hard and soft_QCD dynamics. 

In the Color Singlet Model El (CSM) a colour singlet cc pair is produced 
in a conventional pQCD process. However, this model gives a cross section 
more than one order of magnitude lower than the observed one for high-pj^ 
charmonium at the TevatronH. This surprisingly high rate can be reproduce] 
well by the Color Evaporation Model (CEM)0 and the Soft Color Interaction! 
model which not only include conventional pQCD but also aspects of non- 
pQCD dynamics. In essence, a perturbatively produced cc pair in a color 
octet state can be transformed into a color singlet state through soft gluon 
exchange. 

In CEM a simple statistical color factor 1/9 gives the probability that 
a cc pair is in a singlet state. For pairs with mass below the threshold for 
open charm production, nice < Smu, the fraction giving a J/ip is given bw 
an additional non-perturbative parameter pj/^ =pj 0.5. We have earlier foundH 
that CEM can replace the Pomeron-based modelH to describe photoproduction 
of J/ip and D. As a result of our further investigation of the CEM model, we 
here present preliminary results for the xp and p±^ distributions of J/^ in 
hadroproduction. 

The total cross section for J/^ in hadroproduction was calculated ana- 
lytically using leading order (LO) QCD, where the contributing processes are 
gg — > cc and qq cc. Our results in fig. |l| show an energy dependence in 
agreement with data. However, an overall normalisation factor K = 3.3 was 
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Figure 1: Total cross section of J/?/) and open charm (D) in hadroproduction as obtained 
with CEM based on the LO pQCD subprocesses and including a factor K = 3.3. The curves 
o'open sup and cTopen inf s-^e the Contributions to open charm for nice above and below 2m£) , 
the latter for cc in a color octet state. 

needed in order to describe the data. This may be attributed to either higher 
order pQCD corrections that were not included in our calculation, or to some 
soft dynamics that could also play an important role. 

In order to understand the origin of this JiT-factor, we studied also some 
important differential distributions, namely xp (longitudinal momentum frac- 
tion) and p± of J/tp in fixed target experiments. For this purpose we have 
added the Color Evaporation Model to the PythiaD Monte Carlo program so 
that complete events can be simulated. The LO pQCD processes are the same 
as described above with the charm mass included in the matrix element. The 
importance of higher order pQCD contributions was estimated by the initial 
and final state parton showers where a cc pair can be produced. In this case 
we include all 2 ^ 2 processes (gg — > gg, qg qg, qq — > qq etc) as hard 
scattering basis, except those producing a cc to avoid double counting. 

Our results are shown in fig. eL where the xp and p± distributions of J/i/j 
are compared to fixed target datall. With a charm mass rric = 1.3 GeV, there 
is no need for a K factor, provided the higher order pQCD contributions and 
the soft color effects are included. Thus, both higher order pQCD and soft 
effects are important. 
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Figure 2: Distribjutions in xp and of J/^ in fixed target hadroproduction with Ef,^am = 
800 GeV. Data Q compared to CEM results based on leading order (LO) matrix elements 
and higher orders (HO) estimated by parton showers, and their sum (full histogram). Also 
shown for xp (thin full curve) is the LO result with rric = 1.45 GeV and K = 3.3. 

The higher order pQCD gives important contributions at small xp and 
large p± . The distribution in p± gives further insights on the interplay between 
hard and soft processes. The data can be reasonably well described by a 
gaussian of width 1 GeV, which suggests that both hard and soft contributions 
are present. The model has three basic sources giving tranverse momentum to 
the J/ 4'- The first is from the pQCD processes described by the 2 — > 2 matrix 
elements and by the parton showers. The LO cc process dominates strongly the 



J/V' production at low p±^ , but not at larger p±^ because a larger momentum 
transfer in the 2 — > 2 process leads to Mcc > 2mD and open charm production. 
Instead, the higher order pQCD processes, where typically a gluon from the 
LO 2^2 process splits into a cc pair, is important for the J/ ij) production at 
larger p±^ as can be seen in fig. ^d. 

The two other sources of transverse momentum are due to soft dynamics. 
One is the Fermi motion of partons inside the initial state hadrons, for which we 
use the conventional gaussian distribution with (fc^) = (0.59 Gey/c)^ used in 
Pythia. The other source is the soft momenta that may be associated with the 
soft gluon exchange that neutralize color. We model this source by giving the cc 
pair an additional p±^ in random direction using the same gaussian distribution 
as for the soft Fermi motion. This results in a reasonable description of the 
data as shown in fig. ||b. Thus, both hard pQCD and soft non-pQCD effects 
should be considered in order to understand the observed p±^ distribution. 

In conclusion, we have found that both higher order pQCD and soft non- 
pQCD effects are important for the understanding of charmonium production. 
The soft interactions can be successfully modelled, e.g. by the Color Evapora- 
tion Model which we have made more dynamically explicit by implementing 
it in the Pythia Monte Carlo. Simulated events then give access to more 
detailed information on J /tp production. Our results, which depend on a few 
parameters related to basic unknown features of non-pQCD, provide a reason- 
able description of the overall cross section of J/'0 in hadroproduction, without 
the use of an arbitrary A'- factor. The differential distributions in xp and pj_ 
of the J /-ip are also well reproduced. 
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